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ABSTRACT 



Context. Weak lensing applied to deep optical images of clusters of galaxies provides a powerful tool to reconstruct the distribution 
of the gravitating mass associated to these structures. 

Aims. We use the shear signal extracted by an analysis of deep exposures of a region centered around the galaxy cluster 
ABCG 209, at redshift z ~ 0.2, to derive both a map of the projected mass distribution and an estimate of the total mass within 
a characteristic radius. 

Methods. We use a series of deep archival R-band images from CFHT-12k, covering an area of ~ 0.3 deg 2 . We determine the 
shear of background galaxy images using a new implementation of the modified Kaiser-Squires-Broadhurst KSB+ pipeline for 
shear determination, which we has been tested against the "Shear TEsting Program 1 and 2" simulations. We use mass aperture 
statistics to produce maps of the 2 dimensional density distribution, and parametric fits using both Navarro-Frenk- White and 
singular-isothermal-sphere profiles to constrain the total mass. 

Results. The projected mass distribution shows a pronounced asymmetry, with an elongated structure extending from the SE 
to the NW. This is in general agreement with the optical distribution previously found by other authors. A similar elongation 
was previously detected in the X-ray emission map, and in the distribution of galaxy colours. The circular NFW mass profile 
fit gives a total mass of A/200 = 7.Tt 2 '_7 10 14 Mq inside the virial radius r2oo = 1.8 ± 0.3 Mpc. 

Conclusions. The weak lensing profile reinforces the evidence for an elongated structure of ABCG 209, as previously suggested 
by studies of the galaxy distribution and velocities. 

Key words. Gravitational lensing: weak lensing - Galaxies: clusters - Clusters of Galaxies: individual (Abell 209) 



1. Introduction 

The study of the origin and structure of clusters of 
galaxies occupies a central place in the current efforts to 
understand the origin and evolution of the Universe. In 
particular, the determination of the mass distribution of 
clusters can prove crucial for verifying the existence of 
dark matter in the Universe and eventually to determine 
its abundance and dynamical evolution. 
Cluster masses have traditionally been derived through 
the virial analysis of the velocity dispersion of cluster 
gal axies, with the assumption of dynamical equ ilibrium 
(eg. lGirardi et alJll998t[Girardi fc Mezzettill200lh . and/or 
from the X-ray temperature of t he hot intracluster gas, 
assuming hydrostatic equilibrium (jRosati et al. 2002, and 



references therein). Since using these methods one has 
to assume the dynamical and hydrostatic equilibrium, 
these mass estimates are affected by the ignorance 
about the dynamical state of the cluster. Weak lensing 
analysis offers a unique opportunity to determine the 
cluster mass distribution without such assumptions on 
its equilibrium, as the effect is due to the gravitational 
deflection of the light that is dependent solely on the 
distribution of matter. In particular, one application 
of weak lensing analysis is the detection of the weak 
shear around galaxy clusters, yielding an estimate of 
the total cluster mass and allowing a full mass recon- 
struction of mainly l ow (0.2 ^ z ^ .5) re ds hift clusters 



(Clowe et all l2004t iBardeau et ail 120051 : IClowe et al.1 



2006>J) . As a possible target for a weak lensing analysis 



* This project has been partly supported by a Marie 
Curie Transfer of Knowledge Fellowship of the European 
Community's Sixth Framework Programme, under contract: 
MTKD-CT-002995 COSMOCT. 



the galaxy cluster ABCG 209 is particularly interesting, 
because the photometric and evolutionary properties 
of its g alaxy populations have been already thoroughl y 
studied (jlvlercurio et al.ll2004L E)03aHHaines et al.ll2004j ). 
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In this paper we present the weak lensing mass recon- 
struction of the galaxy cluster ABCG209 at z = 0.21 



(Mc rcurio et alJ I2003U and references therein) using 



archival wide-field i? -band imaging. A BCG 209 is a rich 



(richness class R=3, Abell et al. 19891 ). X -ray luminous 
(L x ( 0.1-2.4keV)~ 14 10 44 erg.s 
11999 1. moderately hot (T x 
and massive cluster ( Dahle et al 



Ebeling 



2002: 



et al 



10 k eV. iRizza et alJll998h 



Mercurio et al 



2003bl) . It is characterized by the presence of substruc- 
tures, which is shown by an elongation and asymmetry 
in the X-ra y emi ssion maps, with two main clumps 
( Rizza et al. 1998). Moreover, the young dynamical 
state is indicated by the presence of a radio halo 



( Giovannini et al. 19991 : Venturi et al. l2006h . which has 



been suggested to be the result of a recent cluster merger, 
through the a cceleration of relativistic particles by the 
merger shocks (|Ferettill2002h . 



The plan of the paper is as follows. In section 2 we will 
shortly review the data reduction procedures. In section 3 
we will present the software pipeline which we have devel- 
oped t o perform this analysis, based on the KSB + algo- 
rithm ( Kaiser et al. 1995 : Luppino fc Kaiser 1997tl . There 



exist many variants of th is algorithm (see iHevmans et al 



2006; iMassev et al. 2006, for a presentation and compari- 



son of some implementations) , so we will describe in some 
detail our particular implementation (the OACt pipeline) . 
In section 4 we will describe the preparation of the galaxy 
catalogue, and the technique adopted to extract the shear 
information. Section 5 will be devoted to present the mass 
estimated and a mass aperture reconstruction obtained 
from the shear maps. In section 6 we will compare the 
mass maps with the galaxy distributions and the X-ray 
maps. Section 7 will present the conclusions. 

2. Data description 

A detailed description of the da ta and reduction t ech- 
niques has been given elsewhere (jHaines et al.l 120041 1 . so 
here we will only summarize the main steps. The data 
were obtained from the Canada-France-Hawaii telescope 
(CFHT) science archive (PI J. -P. Kneib), and comprise a 
wide-field R-band image centered on the cluster. The ob- 
servations were made on 14-16 November 1999, using the 
CFHT12K mosaic camera, an instrument made up of 12 
4096 x 2048 CCDs, set at the prime focus of the 3.6-m 
CFHT. The CCDs have a pixel scale of 0.206", resulting 
in a total field of view of 42 x 28 arcmin 2 , corresponding 
to 8.6 x 5.7 hjQ Mpc 2 at the cluster redshift. The total 
exposure time is 7200 s, made up of twelve 600s exposures, 
jittered to cover the gaps between the CCDs. 
Standard IRAF tools were used to bias-subtract the im- 
ages, using bias exposures and the overscan regions of each 
CCD, before flat-fielding using a superflat made up of all 
science images from the same observing run, registering 
and Co-adding the images. The resultant images have a 
median seeing of 0.73". 

The photometric calibration was performed in the the 



Johnson-Kron-Cousins photometric system, using obser- 
vations of ~300 secondary standard stars (14 <R< 17) in 
fields 6 and 7 of Galadi-Enriquez et al. (2000), resulting 
in a zero-point uncertainty of 0.005 magnitude. 
For the weak lensing analysis we have masked: (i) satu- 
rated stars and hot pixels; (ii) a 1' (ie. 300 pixels) border 
all around the field, where the point-spread function (here- 
after PSF) is too complex to be properly modelled (eg. 
concave and/or varying too rapidly on small scales); (iii) 
CCD gaps (ie. area covered by several CCD on stacked 
images), where the PSF is also too complex. 
The data set does also contain B-band imaging of the field 
with seeing and depth significantly lower than on the R- 
band image. On this B-band image, the bright star den- 
sity is too low to interpolate the PSF, so it is not used to 
compute the shape parameters, but to help distinguish be- 
tween cluster, foreground and background galaxies. Using 
an algorithm which takes into account the R-band mag- 
nitude , B-R colour and the local density, lHaines et al 



(|2004j) attribute to each galaxy a probability of belonging 
to the cluster. Given the lack of spectroscopic information, 
this probability is the most accurate information we have 
to discriminate among cluster and field galaxies. 



3. Weak gravitational shear estimate 

We extract the shear signal from observed polarisation 
of background galaxies, corrected for the effects of the 
PSF, via the standard K SB+ method (jKaiser et al.lll995 : 
Luppino fc Kaiserl Il997j l. improved as described in the 
following sections. We have blind tested our pipeline 
on the Shear-TEsting - Program (STEP) s imulated data 



([Hevmans et al.l l2006t IMassev et al.1 120061 1 . where alto 



gcther 16 different weak lensing pipelines have been tested 
and compared. For reasonable PSFs, as the one sampled 
by the images presented in this paper turned out to be, 
the shear we measure <y meas - is a linear function of the true 
(ie. simulated) one j tTue ; 



yncas. = ( m + 1)^™ + C 



(1) 



where m is the calibration bias and C is a systematic 
effect, mainly describing the anisotropy of PSF residuals. 
Both m and C depend on the PSF, but contain a constant 
factor intrinsic to the pipeline, which can be subtracted. 
After subtraction, one has: 



-0.05 < m < 
-0.02 < ICI < 0.02 



(2) 



For the current data set these systematics are well below 
statistical errors. These statistical errors are mostly due 
to the intrinsic distribution of galaxy polarisations and to 
the background galaxy density. 

We present in sections 13.11 and 13.21 the basic KSB+ for- 
malism, and in section |3~31 the details specific to the OACt 
pipeline. 
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3.1. Polarisation as a shear estimator 

The observed polarisation of a galaxy offers an estimate 
of the local gravitational shear, and can be defined us- 
ing the weig hted quadrupole m oments of the brightness 
distribution (Kais er et al 1 ll995ft : 



J d 2 9W{6)I(6) 6ji 
J d 2 9W(9)I(9) 



(3) 



where I is the surface brightness of the object, 9 is the 
angular distance from the object center and W is a window 
function, whoise introduction is necessary to reduce the 
shot noise to a reasonable level at large distances from the 
centroid. Note that the indexes i and j are symmetric, so 
we have Q12 = Q2i- Using the weighted quadrupoles, we 
define the complex polarisation e = e\ + i ei as: 



1 



Qw + Q22 



Qn — Q22 

2Ql2 



(4) 



For an elliptical object and a constant window function 
W(9) — 1, e is simply related to the axis ratio (3. Defining 
a position angle 9 of the major axis, measured counter- 
clockwise from the x axis, one obtains: 



1 



1 



( cos 29 
/J2 I sin 29 



(5) 



In the weak lensing limit (defined by: I7I <C 1, where 7 = 
7i + ij2 is the complex shear field), and assuming that 
the real (ie. unobserved) polarisation distribution has a 
null average, 7 is directly related to the average observed 
polarisation, 7 w e/2. 

3.2. The KSB method 

The current KSB+ method is the result of a se- 



ries of successive imp rovements (jLuppino fc Kaiserlll997 



Hockstr a et a .1119981) of the original method proposed by 
Kaiser et al.l ( 19951 ). It provides a gravitational shear es- 
timate by first-order subtraction of the PSF smearing 
and shearing from the galaxy polarisation. The 2D vec- 
tor KSB+ shear estimator of a single galaxy 7 is given 
by: 



where we have adopted the standard convention on the 
summation rule of indices, and q is the anisotropic com- 
ponent of the PSF, _P sm is the smear polarisability tensor, 
and P 7 is the pre-seeing shear polarisability tensor, the 
latter being defined as: 



pi — p sti 



PSF 



PSF 



(7) 



Here P sh is the shear polarisability tensor, and the sub- 
script "PSF" signals that the quantity is computed for the 
PSF. In equation^ q and P~* depend on the PSF and are 
estimated from the images of the surrounding stars. 
The actual prescription to estimate q, P 7 and the PSF- 
subscripted tensors, the choice of the window function W 



in equation [21 the algorithm of pixelised summations, and 
finally the approximations, vary from one implementation 
of the method to another. Our algorithm is described in 
the following section. 

3.3. The OACt pipeline 

3.3.1. Smoothing radius, significance, Window 

function, centroid and Summation algorithm 

We determine the significance v and the optimal smooth- 
ing radius r g of each object, v and r g are defined as 
usual in weak lensing: when convolving the image with 
a Mexican-hat filter, the radius r g is the smoothing radius 
for which the object has the best signal-to-noise ratio, and 
the significance v is this best signal-to-noise ratio. 
The window function W of equation [3] is taken to be a 
circular Gaussian centered on the weighted centroid, hav- 
ing a standard deviation equal to r g . The weighted cen- 
troid, computed iteratively, is the point for which weighted 
dipoles are equal to zero: 



d 2 9W{9)I{9) 6i = 



(8) 



Due to the finite size of pixels, all the integrals are re- 
placed by discrete sums with steps of 0.25 pixel in x and 
y directions, truncated at a distance of Ar g . The flux at a 
given position is the linear interpolation of the flux in the 
four nearest pixels, and the background flux in a pixel is 
estimated by SExtractor (more details on the actual pa- 
rameters adopted are given in section l4~Tj) . 



3.3.2. Shape computation 

We compute the ellipticities, smear polarisability tensors 
and shear polarisability tensors of background galaxies as 
described above (section 13. 2p . We do the same for stars, 
except that we do not use the smoothing radii r g of stars 
themselves, but we perform the calculation for every sin- 
gle 0.1 pixel-bin of r g in the range 0.95 < r g < 9.05. This 
bin width of 0.1 pixel is chosen to be much smaller than 
the accuracy of the measured r g (this accuracy is of few 
tenths of pixel), and the range is given by extrcmum val- 



(6) ues of r g in usual data. This gives a total of 81 r g bins. 



After computation of the shape parameters and shear es- 
timator, we discard the objects for which: (i) the standard 
deviation on the centroid is larger than 0.2 pixel; (ii) the 
polarisation |e| is greater than 1; (iii) the shear estimator 
7 1 is greater than 2; (iv) the trace of the P~* tensor is 
lower than 0.2. For the current data set, this corresponds 
to 5% of the background galaxies. Even on noisier data 
sets this usually corresponds to less than 10%. 

3.3.3. Subtracting the smearing and shearing effects of 
the PSF 

The corrective factors of a given galaxy are all computed 
in the r g bin corresponding to the r g of the galaxy itself. 
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In other words, we interpolate PSF properties from stars 
in every single r g bin independently for each galaxy. In 
all the P sm * and P sh * tensors of stars, the off-diagonal 
terms are completely dominated by shot-noise, and typi- 
cally they are more than one order of magnitude smaller 
than the diagonal terms. For this reaso n, following similar 
implementations ( Hoekstra et al.ll2000l ). we systematically 
neglect these off-diagonal terms. Also the non-diagonal 
terms of the tensor P 1 are extremely noisy, and we then 
approximate this tensor as a scalar equal to half its trace. 
These approximations imply that the vector q of equa- 
tion \E\ (ie. the anisotropic component of the PSF) is the 
interpolation of: 



psm* 
aa 



and 



ps 



3sm* aa 



(9) 



(10) 



Finally, we compute the shear estimator following this 
simplified version of equation EJ 



la - pi 



(ii) 



3.3.4. Mapping PSF properties 

From the star catalogues, we interpolate the 4 ratios: 
K/PaT) and (P a \*/P a T) (with a = 1 or 2) over the 
entire region (where the asterisks refer to parameters mea- 
sured on stars). The former two terms give q (equation [5]) 
and the latter two give P 1 through equation [Til] For each 
term, we fit a 2-dimensional, 2-degree polynomial inde- 
pendently on each CCD (and also independently in each 
r g bin). 

As we will see in the following, in order to keep a rela- 
tively high star density, we use a rather permissive star 
selection criterion. As a consequence, our star catalogue 
is contaminated by small galaxies. In order to reject these 
small galaxies when fitting PSF properties, the fits are 
iterative: after the first fit is performed, we reject all ob- 
jects with at least one residual at more than 3tr, and we 
continue to iteratively perform new fits until convergence. 
Typically, this procedure converges after 2 or 3 iterations. 



4. Detection of the weak lensing signal 

To extract the weak shear information from the reduced 
and calibrated CHFT12k images presented in section[2j we 
first build the star-field and the background-galaxy-field 
catalogues. Background galaxies contain the weak lens- 
ing signal, smeared and sheared by the PSF, while stars 
are measures of the local PSF. We first detect all the ob- 
jects within the field and then select those relevant for 
weak lensing. Stars are selected according to their sizes 
and their magnitudes, while the background galaxies are 
selected by cross-checking our relevant object catalogue 



with the galaxy catalogue of lHaines et al.l (|2004l ) . The lat- 
ter contains all galaxies within the field and assigns to each 
a probability of belonging to the field rather than to the 
cluster itself. All these steps are described in section 14.11 
We then compute shape parameters of stars and back- 
ground galaxies (section l3.3.2p . using the stars to map the 
PSF (as described in section l3.3.4p . and finally compute a 
shear estimator for each background galaxy. 



4.1. Star and background galaxy catalogues 

We detect all the objects on the image using SExtractor, 
with very low thresholds. We also get a large proportion 
(~ 50%) of spurious detections but these are rejected 
later. We demand detected objects to have at least 5 
pixels above 1.5 a, where a is the standard deviation of 
the local background (ie. detect_thresh and ANAL- 

YSIS_THRESH set to 1.5, DETECT_MINAREA set to 5), 

where the local background is estimated with keywords 
back_size set to 70 and back_filterize set to 5. 
SExtractor is also used to measure the flux and the 
half-light-radius of each object. We then compute r g 



and is, as described in section 1331 We then successively 
remove from the catalogues : (i) the objects with v < 3; 
(ii) the objects with at least one neighbor nearer than 
3(r g + t 3 (neighbor)); (iii) the objects with at least one 
pixel belonging to the masked area within an aperture 
of 3r g (the masked area is described in section [5]); 
(iv) the objects with r g lower than the local smoothing 
radius of stars. Finally, we compute the shapes of the 
remaining objects and apply shape cuts, as described 
in section 13.3.21 After having performed these steps, we 
have a catalogue of ~ 30000 detections, containing all the 
weak-lensing-relevant background galaxies but also stars, 
cluster galaxies and a large proportion of spurious events 
due to the low thresholds used with SExtractor. 
We build a star catalogue with a loose selection based 
on 5 parameters: the significance v\ the magnitude R\ 
the surface brightness in the central pixel i? max ; the 
half-light-radius r^, and the smoothing radius r g . We 
demand: v > 10; R < 24.0; R + 2.6 < i? max < R + 3.3; 

1.3 < ^4^r < 1.7. This returns 
uniformly dis- 



1.9 < < 2.5; 

1 pixel ' 



1 pixel 

1588 objects (~ 1.5 objects. arcmin -2 ) 
tributcd in the field. We define this catalogue loose 
because v, r g and of stars show variations of ~ 30% 
through the field with the PSF (while we apply these 
constant cuts) . Thus, these cuts result in a star catalogue 
which contains 20 — 45% of non-stellar detections, most of 
them being small galaxies. These fake stars are iteratively 
rejected during the PSF property fits, as described in 
section 13.3.41 At the end of iterations, we are left with 
~ 0.9 stars. arcmin -2 , as shown in figure [3] 
We select background galaxies from the remaining objects 
in two steps. First, we reject fake detections and cluster 
galaxies by cross- checking the remain ing objects with 



the catalogue of lHaines et alj (|2004T) . which includes 



all the galaxies within the field, together with their 
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accurate photometry and the probability for each galaxy 
of belonging to the field rather than to the cluster itself 
(see section [2} . We include in our final catalogue those 



galaxies identified in the catalogue of Haines et al.l (|2004) 
marked as having a probability larger than 80% to belong 
to the field. Second, we reject most of the foreground 
galaxies with the cut R > 21 (since the R-band magnitude 
is the only information we have at this stage). In order to 
reject too faint background galaxies for which the shape 
parameters are not reliable, we cut the ~ 3% faintest 
galaxies for which R > 25.5. 

These cuts are optimised to give a background galaxy 
catalogue with low foreground contamination, while 
keeping almost all the relevant background galaxies. 
They are illustrated in figure [fl where the star sequence 
is clearly visible in red, while background galaxies are 
in blue. The final background galaxy catalogue contains 
16708 galaxies (16.7 galaxies. arcmin -2 ). When combining 
individual galaxy shears to produce shear maps, mass 
reconstructions or density profile fits (see sections 14.31 
and [I]), we weig ht galaxies ac c ordin g to their significance 



as proposed by I Clowe et al. ( 20061 ): we do not consider 



objects with v < 5, while for v > 5 the weight is set equal 
to min(j/ ; 40). In other words, for a given galaxy i, the 
weight Wi is defined as: 

Wi = ; if Vi < 5 

Wi = Vi ;if 5 < Vi <40 (12) 
w t = 40 ; if Vi > 40 

The weighted galaxy number is EjiDj /max(iUj) = 
EjiUj/40 = 8816 (8.8 weighted galaxies. arcmin~ 2 ). 



4.2. Testing the PSF subtraction 

Figure [2] shows the star polarisations over the field before 
and after the anisotropic correction (ie. the subtraction 
of J^i^aTli in equation fTTjl . One can notice the lack of 
large scale correlations due to the PSF anisotropies, after 
the subtraction. The correction also reduces the average 
amplitude and anisotropy of the PSF, as it is clear from 
Figure [31 The corrected distribution of ellipticities is more 
isotropic, and has also a smaller scatter dispersion. 

4.3. Combining the individual shears, building shear 
maps 

To get a shear map, we divide the field into 17 x 11 square 
cells with an overlap of 50% (ie. 50% of the galaxies in 
one cell belong only to this cell, while the remaining 50% 
belong also to at least one of the 8 neighbouring cells) . In 
each cell we average -y according to the weighting scheme 
described in section |4~T1 The resulting shear map is shown 
in Figure |U One can see a characteristic pattern of in- 
creased tangential shear, which coincides with the central 
region of the cluster, as defined by the optical distribu- 
tion of galaxies. This visual impression is confirmed by 
the mass aperture map, as we will see in the following 
sections. 
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Fig. 1. Top (bottom respectively) panel: magnitude - 
central magnitude (optimal smoothing radius r g - magni- 
tude resp.) diagram for all objects detected in the image. 
Objects selected as possible stars are in red, those selected 
as possible background galaxies are in blue. Other objects 
are in black. 



Fig. 2. Star polarisations in the field. Top panel: be- 
fore any correction (ie. measured polarisations). Bottom 
panel: after the anisotropic correction. In both panels, 
the vectors are oriented along the major axis of the 
ellipsoid, their length being proportional to the polar- 
isation: \e\ = V(Qn - Q22) 2 +4Q 2 12 / (Q11 + Q 22 ). The 
scale shown in the upper right corners is |e| = 0.03. The 
straight lines show the (masked) regions covered by several 
unstacked images. 
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Fig. 3. Distribution of stellar polarisations. Upper panel: e\ against e2 of stars before (black) and after (red) the 
anisotropic correction. Bottom panel: |e| distribution (same colours). 
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x (in pixels) 



Fig. 4. Shear map of the whole field obtained with an overlap of 50% (ie. 50% of the field area is accounted into more 
than one cell). The scale is given by the thick line in the upper right corners : 7 = 0.05. 
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5. Mass distribution 

There are quite a few different methods to deduce the 
mass distribution from the individual shears of back- 
ground galaxies (see eg. iBartelmann fe Schneider 2001 , 
for a review). A class of these methods makes use 
of the 2D smoothed shear maps to get the projected 
surface density distribution: ex amples of this class are 
the mass reconstruction meth od ( Schneider fc Seitdll995t 
Seitz fc _S chncidcr| 199a . l200ll ) , or the mass aperture tech- 
nique ( Schneiderl Il996l ). These methods suffer from the 



sheet mass degeneracy: they provide a reliable reconstruc- 
tion of the surface density, but this is known except for 
a constant value. The alternative is to assume a priori 
a given mass profile, so that the degeneracy is removed, 
then attempt a parametric fitting. These latter methods 
are useful to get an estimate of the total mass of the clus- 
ter. Both methods suffer from possible systematic effects, 
for instance the contamination from background sources. 
In section 15.11 we fit two different parametric spherical 
density profiles: a Singular Isot hermal Sphere (SIS) and 



a Navarro- Frenk- White (NFW) ([Navarro et al.l fl996). In 



section [5~2"l we map the mass distribution using mass aper- 
ture statistics. 

In order to derive actual values of the mass from the 
shear parametric, we have to know the critical surface den- 
sity: 

(D ls /D s 



S C rit — 



(13) 



where D\ s , D\ and D s are the lens-source, observer-lens 
and observer-source distances respectively. We do not have 
however any redshifts of the background galaxies, so we 
approximate their redshift distribution by a single me- 
dian value (ie. we adopt a source plane approximation). 
As described in section I4.1( we restrict the background 
galaxies to the range 21 < R < 25.5: in this magnitude 
range, we ca n assume the median redshift to be z ~ 1 



(jGavazzi et alj |2004) . Then the critical density at the red- 
shift of ABCG209 (z = 0.21) becomes: 



: crit ~ 1.62 x I0 14 M Q .arcmhr 
~ 3.91 x 10 15 M Q .Mpc~ 2 



(14) 



for the basic ACDM c osmological model derived from 
the 3-year WMAP data (jSpergel et al.ll2006l ): fl M = 0.27, 
Oa = 0.73, H = 70km.sec" 1 .Mpc- 1 and w = -1. 



5.1. Parametric model fitting 

For both profiles we have performed a \ 2 minimization of 
the tangential shear: 



A' 



53 Wi ~ 7modcl 



(15) 



where -f T is the tangential shear (ie. the shear projected 
along the direction orthogonal to the line connecting the 



galaxy position to the cluster center), Wi is the weight (as 
defined in equation[12]) and Ti is the distance to the cluster 
center of mass. The latter is defined as the point for which 
the signal-to-noise of the -f T profile attains its maximum 
value: 

o N 

^(a,5)=^2w i7 f(a,8) (16) 

i=l 

By maximising eq. 1161 we find the center at 
a = Ih31m51.6s ; S = -13°36', offset by about 36" 
from the cD galaxy, as shown on figures [S] 
For the NFW fit we have either assumed that each galaxy 
provides an estimate of the underlying shear field, (so the 
summation in equation [T5] extends to the full sample) or 
we have binned the data and minimised over the averages 
in the bins. We have verified that both procedures give 
compatible results. 

We exclude from the fits galaxies lying within an inner 
central region of the cluster, where the weak lensing 
approximation (n < 1) is not everywhere valid, and 
outside an outer region, where the statistics are very 
weak. 

The actual fitted region lies within the following bounds: 
V < 9 < 10' (0.2 Mpc < R < 2Mpc) around the cluster 
center of mass. 

5.1.1. Fitting a NFW profile 

The NFW profile has often been used as a good fit o f 
numerically simulated halos (jNavarro et al.l Il995l 119971 ) . 
Although this fit was originally made only for simulated 
halos in standard CDM models, it turned out to be a good 
fit for halos which formed in ACDM models. The mass 
density of the NFW profile is described by: 

ScPc 



p(r) 



where: 5 C 



(r/r s )(l 
200 



r/r 



and: 



Pc = 



3 ln(l + c) - c/(l + c) 
3H 2 (z) 



BttG 



(17) 



(18) 



(19) 



In the equations above r is the distance to the cluster 
center, r s is the scale radius, H(z) the Hubble parameter 
at the redshift of the cluster and c = r2oo / r s is the con- 
centration parameter, where r2oo is the virial radius. The 
NFW density profile is shallower than the SIS profile near 
the center but steeper in the outer parts. The total mass 
inside the radius R, shown on figure is: 



M(< R) = 4tt S c p c r 3 s ln(l + R/r s ) - 1 + 



1 



1 + R/r s 



(20) 

Exact expressions for th e tangential shear due to th is mass 
distribution are given by lWright fc Brainerd ( 2000l ). It has 
two highly correlated degrees of freedom: r s and c. From 
a non-linear least-squares fit (Levenberg-Marquardt) we 
obtain: 

O ,1+3.1 

(21) 
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This corresponds to virial radius r 2 QQ and mass M\ 



200- 



7"200 

M 200 



1.8lt°J°Mpc 
7.7±Hl0 14 M n 



(22) 



5.1.2. Fitting a SIS profile 



For a SIS model, the density and shear profiles are respec- 
tively given by: 



p(r) = 



2irGr 2 



c 2 D s 29 



(23) 
(24) 



where a v is the velocity dispersion, the only free parameter 
in this model. The total mass inside the radius R (shown 
on figure [HJ is: 



M(< R) = -<J 2 V R 
A linear least-squares fit provides: 

<t„ = 924±84km.s _:l 



(25) 



(26) 

The previous analysis from Dahle et al. ( 20021 ) gives a sig- 
nificantly lower value (a v = 680^^3° km. s _1 ), but the au- 
thors adopt a less restrictive criterion to eliminate cluster 
galaxies before the weak lensing analysis, so their weak 
lensing signal, and consequently a v , could be underesti- 
mated. Although the total mass of a SIS model diverges, 
we can however show the integrated mass inside a given ra- 
dius R (Figure^). At the virial radius R = r 200 ~ 1.8 Mpc, 
as provided by the NFW profile fit (previous section), the 
mass inside the radius is almost the same both for SIS and 
NFW profiles: M(< R) = M 20a ~ 7.7 x 10 14 M Q . 
In Figure[5]we show the total mass inside the radius R, ac- 
cording to the two different fits we find. Note that the ±lcr 
region of NFW is much larger than in the SIS case. This 
is a direct consequence of equations [20] and [25] in the SIS 
model M (< R) is proportional to c 2 and so the confidence 
region has a width of ~ ±20%, while in the NFW model 
M(< R) depends on S c r^f(r s ,R). This makes the ±ler re- 
gion to have a minimum ~ ±50% width for R ~ 1.3 Mpc. 



5.2. Mass aperture statistics 

The mass within an aperture radius r can be obtained 
directly from the shear using the aperture de nsitometry 
statistics £, defined as in iFahlman et"aT] (|l994 : 



C(Ri,R 2 ) = fi(Ri)-H(Ri,R 2 ) 



l-R\lRlJ Rx 



fB.2 

{j T ) r d lnr 



(27) 
(28) 



The mass inside the annulus defined by (R\,R 2 ) is con- 
nected to the above quantity by: Af ap = 7rr 2 ££ cr it- This is 
in fact a lower limit on the true mass, but only dependent 
on the tangential shear. Thus, it is not affected by residual 
B-modes and gives a non-parametric representation of the 
matter density at a given point. 



Following the approach of Schneider ( 19961 ). we build 
the (2 D) mass aperture statistics by computing the aper- 
ture densitometry at each point of a grid on the field. The 
mass density M ap at a given point is given by: 



Map = 



Sj Wi Qi 



(29) 



where the sum extends over all the background galaxies 
and Q is a smoothing-weighting function. We use the win- 
dow function Q that maximises the t otal signal-to-noise 
ratio (S/N) for a SIS (|Schneiderill9"96h : 



Q{x) 



6nx 2 x (1 




;x < 
;x > 



(30) 



where: x = r/r an , r being the distance to the center of the 
annulus, and r an = 5 arcmin (equivalent to = 1500 pixels 
= 1 Mpc) is the external radius of the annulus. The only 
free parameter, ie. the filter scale r an , is chosen as a com- 
promise between different needs. On one side, r an needs 
to be as small as possible to avoid to wash out small scale 
structures. On the other hand the annulus needs to be 
large enough to encompass a significant number of back- 
ground galaxies (typically 1000) everywhere on the field, 
in order to get a good and spatially stable S/N. Figure [5] 
shows the mass aperture density map of the ~ 100 arcmin 2 
patch centered on ABCG 209, together with the cluster 
center of mass defined by equation llGI The mass aperture 
statistics show the elongated structure of the mass dis- 
tribution. The mass distribution is far from the circular 
symmetry as assumed by SIS and NFW profile fits. This 
explains why these fits give only a rough estimate of the 
total mass, as discussed in section [7] 

6. Comparisons with other observations 

6.1. Comparison with the galaxy distribution 

Figure E] shows the galaxy distribution (shown as the 
grayscale-filled black contours) as represented by the sur- 
face density of R < 23.0 cluster galaxies (ie. after correct- 



ing fo r background/foreground galaxies, see lHaines et al 



2004). The galaxy distribution is strongly elongated in the 
same SE-NW axis as observed for the weak-lensing recon- 
structed projected mass distribution. The center of the 
galaxy distribution appears somewhat offset to the NW 
by about 1' with respect to the center of mass obtained 
from the weak lensing reconstruction. There is also a sub- 
structure about 5' to the North of the cluster center, which 
could be connected with the observed substructure in the 
reconstructed mass distribution that extends northwards 
from the central mass concentration. 

The interna l dynam ics of the cluster was studied by 
Mercurio et al.l (|2003af ) through a spectroscopic survey of 
112 cluster members. A high value of the line-of-sight 
velocity dispersion was found, with o v — 1250^g4kms -1 
after removing seven interlopers. Assuming dynamical 
equilibrium, this value of er„ leads to a virial ra- 
dius of R V i r ~ 3.28 ± 0.55 h^ Mpc and a virial mass 
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Fig. 5. Total mass inside the radius R for the two models SIS (in blue) and NFW (in red). Dashed lines are the 66% 



confidence level regions. The black vertical lines show the virial radius i?2oo = l-81_Q'2gMpc given by the NFW fit: 
straight line for the best fit and dashed lines for the 66 % confidence level regio n. The black cross (in the upper right 
corner) shows the measure of M(< R V i r ) and R V i r bv iMercurio et aL ( 2003a ) from internal dynamics, discussed in 
section I6.ll 



of M(< R mr ) = 3.02±g:||, x 10 15 /itoM q in a ACDM 
model, with f2 m =0.27 and f2A=0.73. We report this value 
in Figure O for a direct comparison with the results of 
the NFW and SIS analyses showed in this paper. The dif- 
ference in the estimates of the virial mass and the virial 
radius obtained by the kinematics and weak lensing could 
be due to the presence of substructures which results in an 
over estimate of the velocity dispersion. Another source of 
uncertainties in the kinematics could be the anisotropy pa- 
rameter of the cluster velocity distribution. On the other 
hand, the lensing values could be biased because of the 
elongation of the cluster mass distribution (see section [7]). 

Evidence in favour of the cluster undergoing a dynam- 
ical evolution is found in the form of a velocity gradient 
acting along a SE-NW axis, which is the same preferen- 
tial direction found from the elongation in the spatial dis- 
tribution of galaxies, as well as that of the cD galaxy. 
There is also significant deviation of the velocity distri- 
bution from a Gaussian, with evidence for two secondary 
clumps at z = 0.199 and z — 0.215, which appear spatially 
segregated from the main cluster. These all indicate that 
ABCG 209 is undergoing strong dynamic evolution with 
the merging of two or more sub-clumps along the SE-NW 
direction. 



6.2. Comparison to X-ray emission 

The X-ray data are taken from th e XMM science arc hive 
(Prop #8423, PI. J.-P. Kneib, see lMartv et"aH (|2003l ) for 



an analysis). The observations were made in Jan 2001, 
with an exposure time of 20ksec. The EPIC MOS1, MOS2 
and pn images were combined over the temperature range 
0.5 — 12keV and the resulting spatial distribution of the 
X-ray emission is shown in Fig. [7| by the grayscale-filled 
contours. The X-ray emission is centered on the cD galaxy 
(a =lh31m52.5s, S = -13°36'40", z = 0.2097), making it 
slightly offset (36") from the center of mass determined 
from the weak lensing by equation [TCI The X-ray emission 
is elongated along the same SE-NW direction as seen 
for the weak lensing reconstructed mass distribution, 
the emission being most extended towards the NW. 
There is no evidence of excess X-ray emission from the 
substructure seen in the weak-lensing reconstruction ~ 3' 
to the North of the cluster center. 

From an analysis of a lOksec Chandr a ACIS-I (0.3 



lOkeV) X-ray observation of the cluster, IMercurio et al 
(|2003ah obtained a best-fitting temperature of 
Tx = 10.2+};2keV which, assuming f3 spec — 1, would 
correspond to a v ~ 1300km.s _1 , and is consistent with 
the high value of Ly (0.1 - 2.4keV)~ 14 lO^erg.s" 1 
(jEbeling et al.lll996l ) . This value of the velocity dispersion 
produces a virial mass estimate of 3.3 x 10 15 M Q 
The mass estimates obtained through our weak lensing 
analysis are lower than those based on the X-ray temper- 
ature and galaxy velocity dispersions. In a weak lensing 
analys i s of 3 5 X-ray luminous clusters at 0.15 < z < 0.30, 
Dahle ( 2006) finds a large scatter in the relation between 
the weak lensing mass estimates and the X-ray luminosity, 
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Fig. 6. Comparison of the weak lensing mass reconstruction with the galaxy distribution for ABCG 209. The 
black contours represent the isodensity contours from the M ap statistics, corresponding to mass densities 5 x 
10 13 Mq arcmin -2 (thin dot-dashed curve) to 1 x 10 14 Mq arcmin -2 (tick solid curve), with each contour separated 
by 1 x 10 13 Mq arcmin -2 . The cluster center of mass as defined by Equation 14 is indicated by the red diamond. The 
grayscale filled black contours represent the surface density of R < 23.0 galaxies (background-corrected), corresponding 
respectively to (dashed curve) 2 (solid), 3.3, 5, 7.5, 10.0 and 12.5 cluster galaxies. arcmin -2 . 



producing a mass uncertainty of um = 0.44 dex. In a 
weak lensing st udy of 24 X-ra y lumi nous clusters at 



0.05 < z < 0.31, ICvpriano et all (|2004h found that on 



average the mass estimates based on X-ray temperatures 
and velocity dispersions were 13-27% higher than those 
from the weak lensing analysis. In particular they found 
the discrepancy to be much greater for the most massive 
clusters (Tx > 8keV or a v > 1122 km. s -1 ), where the 
mass excess from the X-ray temperatues or velocity dis- 
persions were 40-75%. They found that the discrepancy 
was largest for the two clusters with the largest X-ray 
temperatures (Tx ~ 13keV) and velocity dispersions, 
which were known to be undergoing a merging event, and 
are far from equilibrium. The high X-ray temperatures 
would then be probably due to recent shocks, and the 
high velocity dispersions due to substructures and the 
complex dynamical situation. 



6.3. Comparison with Galaxy Colours 

The star-formation history of galaxies is known to cor- 
relate strongly with their local environment. In Figure [5] 



we compare the mass distribution with the mean B — R 
colour of the R < 21 cluster galaxy population as a func- 
tion of spatial position. Each galaxy is weighted according 
to the probability that it belongs to the cluster, and then 
the mean colour of cluster galaxies calculated as a function 
of spatial position using an adap tive kernel m e thod. This 
analysis is described in details in lHaines et al.l (|2004t) : the 
mean galaxy colours (and hence their star-formation his- 
tories) are strongly correlated with the dynamical state of 
the cluster, with an alignment of the colours with the main 
SE-NW axis. The region with the reddest mean galaxy 
colours, and hence the oldest stellar populations, is found 
at the cluster center of mass. There are also two regions of 
red galaxies outside the cluster core that are aligned with 
the dark matter distribution, confirming that galaxy evo- 
lution is strongly dependent on the hierarchical build up 
of clusters through mergers. Given the uncertain effects 
cluster dynamics have on the X-ray emission and galaxy 
velocity distributions, maps of the mass distribution based 
on weak lensing analyses provide an important tool for un- 
derstanding the relation between gala xy evolution and the 



underlying dark matter distribution ( Gray et al. 20041 ). 
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Fig. 7. Comparison of the weak lensing mass reconstruction with the X-ray emission for ABCG 209. The red contours 
represent the isodensity contours from the M ap statistics, as in Fig. \E\ The black grayscale-filled contours represent 
the X-ray emission based on the XMM imaging. The contours are logarithmically spaced, adjacent contours indicating 
a factor two change in flux density. 



7. Conclusions 

We have performed a weak lensing analysis of the galaxy 
cluster ABCG 209 through a new implementation of the 
KSB+ algorithm (the OACt pipeline), and we have also 
performed a mass reconstruction using Mass Aperture and 
parametric statistics. We clearly find a measurable weak 
lensing signal, and the comparison with optical and X- 
ray data for this cluster brings some interesting conlu- 
sions. First, the centers of the X-ray emission, dark mat- 
ter, and galaxy distribution all appear offset from one an- 
other, with the center of mass found from the weak lensing 
analysis lying between that of the X-ray and galaxy dis- 
tributions, and all the three centers of mass aligned on the 
main SE-NW axis of the cluster. Such an effect is seen for 



the more extreme Bullet cluster ( Clowe et al. 2004), and 
appears to reflect the different responses of the gas and 
dark matter components to the merger, briging a further 
hint at a cluster merging scenario for ABCG 209. Second, 
we confirm that ABCG 209 is a massive cluster, although 
the mass estimate d by weak lensin g is lowe r than the esti- 
mates obtained by lMercurio et ah (2003a) from the anal- 
ysis of the dynamical properties of the galactic population 
(assuming the dynamical equilibrium). On the weak lens- 
ing side, there are two sources of error not taken into ac- 
count in this analysis, that could explain this discrepancy. 
First, the 2D-mass distribution of the cluster is not circu- 
lar and this fit of a circular profile is possibly not accurate. 
Second, we have not taken into account the uncertainty on 
the critical surface density (equations [Till and [T4"|) . which 
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Fig. 8. Comparison of the weak lensing mass reconstruction with the mean colours of galaxies across ABCG 209. The 
red contours represent the isodensity contours from the M ap statistics, as in Fig. ® The coloured contours indicate 
the mean B — R colour of the R < 21 cluster galaxy population (after statistically correcting for field contamination) 
as a function of spatial position. 



has been computed according to the single-source plane 
approximation at z = 1. Considering these uncertainties, 
the agreement we find among the different mass estimates 
should be regarded as satisfactory. 
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